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Abstract The systematic magnification of background sources by the weak gravitational-lensing effects of foreground matter, 
also called cosmic magnification, is becoming an efficient tool both for measuring cosmological parameters and for exploring 
the distribution of galaxies relative to the dark matter. We extend here the formalism of magnification statistics by estimating 
the contribution of second-order terms in the Taylor expansion of the magnification and show that the effect of these terms was 
previously underestimated. We test our analytical predictions against numerical simulations and demonstrate that including 
second-order terms allows the accuracy of magnification-related statistics to be substantially improved. We also show, however, 
that both numerical and analytical estimates can provide only lower bounds to real correlation functions, even in the weak lens- 
ing regime. We propose to use count-in-cells estimators rather than correlation functions for measuring cosmic magnification 
since they can more easily be related to correlations measured in numerical simulations. 
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1. Introduction 

Gravitational lensing by large-scale structures magnifies 
sources and distorts their images. The systematic distortion of 
faint background galaxies near matter overdensities, the cos- 

• mic shear, has been measured by several groups in the past few 
years (Bacon et al. 2000, 2002; Hammerle et al. 2002; Hoekstra 

■ et al. 2002; Kaiser et al. 2000; Maoli et al. 2001; Refregier et 
al. 2002; Rhodes et al. 2001; Van Waerbeke et al. 2000, 2001, 
2002; Wittman et al. 2000). It was found to be in remarkable 
agreement with theoretical predictions based on the Cold Dark 
Matter model, and has already provided new constraints on cos- 
mological parameters (Van Waerbeke et al. 2001). 

In a similar way, systematic magnifications of background 
sources near foreground matter overdensities, the cosmic mag- 
nification, can be measured and can provide largely inde- 
pendent constraints on cosmological parameters (Menard & 
Bartelmann 2002, Menard et al. 2002). Gravitational magni- 
fication has two effects: first, the flux received from distant 
sources is increased, and the solid angle in which they appear 
is stretched, thus their density is diluted. The net result of these 
competing effects depends on how the loss of sources due to di- 
lution is balanced by the gain of sources due to flux magnifica- 
tion. Sources with flat luminosity functions, like faint galaxies, 
are depleted by cosmic magnification, while the number den- 
sity of sources with steep luminosity functions, like quasars, 
is increased. Thus, cosmic magnification gives rise to appar- 
ent angular cross-correlations between background sources and 



foreground matter overdensities which are physically com- 
pletely uncorrelated with the sources. These overdensities can 
be traced by using the distribution of foreground galaxies. 

Numerous studies have confirmed the existence of quasar- 
galaxy correlations on angular scales ranging from one arc 
minute to about one degree, as expected from cosmic lens- 
ing (for a review, see Bartelmann & Schneider 2001; also 
Guimaraes et al. 2001). However, the measured amplitude of 
these correlations has been systematically higher than usually 
theoretically predicted. 

While cosmic shear can directly be related to observable 
quantities like image ellipticities, the theoretical interpretation 
of cosmic magnification involves several approximations: 

- the luminosity function of the sources is described by a 
power-law over the range probed by the flux limit of the 
observation; and 

- the magnification is assumed to fall into the weak lensing 
regime, i.e. to deviate weakly from unity. Thus, the mag- 
nification can with sufficient accuracy be approximated by 
its first-order Taylor expansion and its deviation from unity 
becomes proportional to the lensing convergence alone. 

While the first assumption is comfortably satisfied, in par- 
ticular for quasars, the validity of the second needs to be veri- 
fied. This is the goal of the present paper. 

Our paper is structured as follows: first, we introduce the 
formalism of the effective magnification and its Taylor expan- 
sion in Sect. 0. We then describe a number of statistics related 
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to the lensing convergence, and evaluate the amplitude of the 
second-order terms which appear in the Taylor expansion. In 
Sect. [3[ we describe the numerical simulations we use to test 
our analytical results and estimate the accuracy of several ap- 
proximations for the magnification. As an application, we in- 
vestigate second-order effects on quasar-galaxy correlations in 
Sect. and we summarise our results in Sect. |J 

2. Formalism 

2. 1. Expanding the magnification 

Cosmic magnification can be measured statistically through 
characteristic changes in the number density of the background 
sources. Along a given line-of-sight, this effect depends on two 
quantities: 

- the magnification factor \i, which describes whether 
sources are magnified or demagnified, depending on 
whether the matter along their lines-of-sight is preferen- 
tially over- or underdense compared to the mean, 

- and the logarithmic slope a of the source counts as a 
function of flux, which quantifies the amplitude of source 
number-count modifications due to flux magnification. As 
mentioned in the introduction, magnification by gravita- 
tional lensing not only increases the observed flux, but also 
stretches the sky, thus the number density of sources on a 
magnified patch of the sky is reduced. The net magnifica- 
tion effect, called magnification bias, depends on the bal- 
ance between the number of sources lost by dilution and 
gained by flux magnification. The steeper the number-count 
function of the sources is, the more pronounced is the mag- 
nification bias. 

If the number-count function of the background sources can 
be described as a power law in a sufficiently wide range around 
the flux limit of the observation, the magnification bias is quan- 
tified by the effective magnification /i Q_1 . It directly expresses 
the changes of the background source density caused by lens- 
ing through the relation 



i(> S,0) = ^- 1 {6) m {> S) 



(1) 



where no(> S) is the intrinsic number-count function of 
sources whose observed flux exceeds S in the absence of lens- 
ing, and n(> S) is the corresponding number-count function 
in presence of lensing. 

The local properties of the gravitational lens mapping are 
characterised by the convergence k, which is proportional to 
the surface mass density projected along the line-of-sight, and 
the shear 7, which is a two-component quantity and describes 
the gravitational tidal field of the lensing mass distribution. The 
effective magnification is related to k and 7 through 

/*— 1 = [(i - «) a - l-yl 8 ] 1- ™ , (2) 

where I7I = (7^ + 7|) 1 ^ 2 is taken as the absolute value of 
the shear. In the weak-lensing regime, both k and \ j\ are small 
compared to unity, and the previous expression can be ex- 
panded in a Taylor series: 

/i"" 1 = l+(a-l) [2k + {2a - 1)k 2 + | 7 | 2 ]+0(k 3 , | 7 | 3 ) .(3) 



Previous studies using analytical formulae for magnification 
statistics focused only on the first-order term of this expansion, 
i.e. they used the approximation /i" -1 « 1 + 2(a— 1) K, which 
potentially causes the amplitude of the effect to be underesti- 
mated. In this section, we investigate the second-order terms in 
the expansion and estimate their contribution. 

In doing so, we first note that n 2 {0) and |7| 2 (0) share the 
same statistical properties (e.g. Blandford et al. 1991), because 
both k and 7 are linear combinations of second-order deriva- 
tives of the lensing potential. The identity of their statistics is 
most easily seen in Fourier space. Since we will only deal with 
ensemble averages of the magnification later on, k 2 and | 7 | 2 
can be combined into a single variable, which we denote by n 
for simplicity. Thus, we can write for our purposes, 



if 1 ' 1 = 1 + 2 (a - 1) [« + an 2 ] + 0{k 3 ) . 



(4) 



Observable effects are due to departures from the mean value of 
the magnification. Therefore, the relevant quantity to correlate 
is = /i Q_1 — Then, up to second order in n 2 , 

the autocorrelation function of the effective magnification is 



4(«-l) 2 [<«(0)«(0 + 0)> 
2a(K(cj>)K 2 (ct> + e))] , (5) 



and the corresponding power spectrum can be expanded in a 
similar way, 



P^-i(s) = 4(a - l) 2 [P K {s) + 2aP fi , 2 (s)} 



(6) 



the power spectrum P^^{s) will be defined in Eq. ( |l6| ) below. 
The last two equations show that the importance of the second- 
order terms in the expansion ^ increases as the number- 
count function of the background sources steepens, i.e. as a 
increases. In the following, we will use a = 2 unless stated oth- 
erwise. This value applies, for instance, to the number counts 
of bright quasars with tub < 19.5 (Pei 1995). For simplicity, 
we abbreviate (SfiSfi) by (////). 

2.2. Second and Third-Order Correlations 

We will now estimate several ^-related statistical quantities 
needed in the Taylor expansion of the magnification. For this 
purpose, we first introduce the k projector such that 



k(0) 



dwp K (w)S[0f K (w),w] 



(7) 



can be written as a weighted line-of-sight projection of the den- 
sity contrast 6 from the observer to the Hubble distance wr. 
The projector is 



WH dw' , f K (u>)f K (w'-w) 
n s {w ) - 



a(w) 



fHW) 



(8) 



where w is the radial coordinate distance, fx (w) is the comov- 
ing angular-diameter distance, ns(w) is the normalised dis- 
tance distribution of the sources, and a(w) is the cosmological 
scale factor. Using Limber's equation, we can then relate the 
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autocorrelation function of k to the dark-matter power spec- 
trum P§, 



(ft(0)ft(</> + 0)) = / dw 



sds 



2tt 



Pi 



fk(w) 



(9) 



where Jo is the zeroth-order Bessel function, and the power 
spectrum P K corresponding to this correlation function is 



A pK w ) n 



fk(w) 



(10) 



As indicated by Eq. (Q), the estimation of second-order 
terms requires the computation of the cross-correlation be- 
tween k and K 2 . We do this by first introducing a three-point 
correlation function for k and then identifying two of its three 
points. As usual, we define the three-point function by 

z K {0i,0 2 ) = (ft(0)ft(0 + 0!)ft(0 + 2 )) . (11) 

Using the k projector defined in (§), we can then write 

z K (9i,0 2 ) = J dwip K (wi) J dw 2 p K {w 2 ) J dw 3 p K (w 3 ) 

x (5[/jc(u;i)0,u>i] 6{f K (w 2 )(4> + 0%), w 2 ] 

x 5[f K (w 3 )(cf> + e 2 ),W3]) . (12) 

Next, we employ the approximation underlying Limber's equa- 
tion, which asserts that the coherence length of the density fluc- 
tuation field is much smaller than the scales on which the pro- 
jector p K varies appreciably. Finally, we insert the expression 
for the bispectrum of the dark-matter fluctuations detailed in 
Appendix and find 



«2_ 

(2tt) 2 " J (2tt) 2 

x B s (ki,k 2 , -ki - fc 2 , w) e ik2 ° 2 t*M , (13) 



where Bg(ki, k 2 , k 3 ) is defined by 

(5(fci)5(fc 2 )5(fc 3 )) = M*i + fc 2 + fcs) B s (k!,k 2 ,k 3 ) .(14) 

Then, using Eq. (|l2|) and identifying two points of the three- 
point correlation function B\ — * 6 2 = 9, we find 



(ft(0)ft 2 (0 + 0)> 



dwp^.(w) 



d 2 fci 
(2^ 



iki-0 f K (-b 



B s (ki,k 2 , -fei - fe 2 , (.15) 



The term (k(<P)k 2 (4> + 6)) is a function of only. Its contri- 
bution P fJ ,. 2 (s) to the power spectrum of the magnification is 
given by the inverse Fourier transform of Eq. (O), which reads 



„ . , f, p 3 Jw) f d 2 t 



s — s — s 



fk(w) fk(w) fk(w) 



(16) 



2.3. Results and predictions 

We can now numerically evaluate the first two contributions to 
the Taylor expansion of the magnification autocorrelation func- 
tion defined in Eq. (|^). As mentioned before, we use a = 2 
here. 



4.5x10 ^ - 



Without Smoothing 



3.0x10 - 




Figure 1. The upper panel shows the amplitude of the two first 
terms of the Taylor expansion of the magnification autocorrela- 
tion, namely (k(4>)k(4> + 6)) (dashed line) and 4(ft(</>)ft 2 (</> + 
0)) (dotted line), using a source redshift of unity. The sum of 
these two terms is shown as the solid line. The lower panel 
details the relative contribution of the second-order term 
for different source redshifts. The figure shows that the lowest- 
order approximation /i k 1+2k misses a substantial part of the 
amplitude of the magnification autocorrelation function. Given 
the accuracy of the bispectrum fitting formula, is accurate 
to -2%. 



For evaluating the correlation functions, we use a CDM 
power spectrum in a spatially flat Universe parameterised with 
Q.0 = 0.3, cr 8 = 0.9, h = 0.7 and T = 0.21. The non-linear 
evolution of the power spectrum and the bispectrum are com- 
puted according to the formalisms developed by Peacock & 
Dodds (1996) and Scoccimarro et al. (2000), see Appendix^. 
The upper panel of Fig. [I] shows the first- and second-order con- 
tributions (dashed and dotted lines, respectively) to the Taylor 
expansion of the magnification for a fixed source redshift of 
z s = 1. The sum of the two contributions is shown by the solid 
line. The figure shows that the contribution of the second-order 
term reaches an amplitude of more than 30% of the first-order 
term on angular scales smaller than one arc minute. According 
to Eq. (S) which describes the Taylor expansion of the magni- 
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fication autocorrelation, we define the relative contribution of 
the second-order compared to the first-order term as 



Rpn{Q) - 



2a( K (c/>)K 2 (0 + fl)) 

<K(»K(> + 0)> 



(17) 



The lower panel shows this ratio in per cent for different source 
redshifts as a function of angular scale. From the lower to the 
upper curves, the source redshifts are 1, 1.5, 2 and 3. For each 
source redshift, the contribution of the second term exhibits a 
similar dependence on angular scale: 

- on scales larger than a few degrees, the contribution drops 
to negligible values, 

- effects become relevant on smaller scales, with a fairly con- 
stant amplitude from a few degrees down to around 10 arc 
minutes, 

- on yet smaller scales, the second-order contribution in- 
creases steeply, due to the non-linear evolution of the den- 
sity field. For sources at redshift 2, the amplitude of the 
second term reaches half of the amplitude of the first term 
below one arc minute. 

Thus, given the amplitude of R^, the correcting term in- 
troduced in Eq. (^) is relevant and must be taken into account 
for describing the magnification autocorrelation with an accu- 
racy better than 30% — 50% on scales smaller than a few de- 
grees. 

So far, we have only investigated the amplitude contributed 
by the second-order term. In order to estimate the remaining 
contributions of all missing terms of the magnification expan- 
sion, we will now use numerical simulations allowing a direct 
computation of /j, as a function of the convergence k and the 
shear 7. 

3. Magnification statistics from numerical 
simulations 

3. 1. The ray-tracing simulation 

For testing the theoretical predictions we performed ray-tracing 
experiments in a Very Large A-body Simulation (VLS) re- 
cently carried out by the Virgo Consortium (Jenkins et al. 2001, 
and see also Yoshida et al. 2001 for simulation details)]] 
The simulation was performed using a parallel P 3 M code 
(MacFarland et al. 1998) with a force softening length of 
^soft ~ 30/i -1 kpc. The simulation employed 512 3 CDM par- 
ticles in a cubic box of 479 /i _1 Mpc on a side. It uses a flat 
cosmological model with a matter density flo = 0.3, a cos- 
mological constant tt\ = 0.7, and a Hubble constant h = 
0.7. The initial matter power spectrum was computed using 
CMBFAST (Seljak & Zaldarriaga 1996) assuming a baryonic 
matter density of 57b = 0.04. The particle mass (m par t = 
6.86 x 10 10 hr 1 Mq) of the simulation is sufficiently small 
to guarantee practically no discreteness effects on dark-matter 
clustering on scales down to the softening length in the redshift 



1 The ray-tracing data are available from T. Hamana on request, 
hamanatk@cc.nao.ac.jp 



range of interest for our purposes (Hamana, Yoshida & Suto 
2002). 

The multiple-lens plane ray-tracing algorithm we used is 
detailed in Hamana & Mellier (2001; see also Bartelmann & 
Schneider 1992 and Jain, Seljak & White 2000 for the theoret- 
ical basics); we thus describe only aspects specific to the VLS 
TV-body data in the following. In order to generate the density 
field between z = and z ~ 3, we use a stack of ten snapshot 
outputs from two runs of the A-body simulation, which differ 
only in the realisation of the initial fluctuation field. Each cubic 
box is divided into 4 sub-boxes of 479 2 x 119.75/i~ 3 Mpc 3 with 
the shorter box side being aligned with the line-of-sight direc- 
tion. The A-body particles in each sub-box are projected onto 
the plane perpendicular to the shorter box side and thus to the 
line-of-sight direction. In this way, the particle distribution be- 
tween the observer and z ~ 3 is projected onto 38 lens planes 
separated by 119.75 /i _1 Mpc. Note that in order to minimise 
the difference in redshift between a lens plane and an output of 
A-body data, only one half of the outputs (i.e. two sub-boxes) 
at z — are used. 

The particle distribution on each plane is converted into the 
surface density field on either a 1024 2 or 2048 2 regular grid 
using the triangular shaped cloud (TSC) assignment scheme 
(Hockney & Eastwood 1988). The two grid sizes are adopted 
for the following reasons: 

- the 1024 2 grid is chosen to maintain the resolution provided 
by the TV-body simulation and removing at the same time 
the shot noise due to discreteness in the A-body simulation. 
Its computation follows the procedure described in Hamana 
& Mellier (2001) and Jain et al. (2000). The corresponding 
outputs will be labelled with large-scale smoothing in the 
following. 

- the 2048 2 grid is also chosen to examine effects of small- 
scale nonlinear structures which are smoothed in the large- 
scale smoothing simulation. We should, however, note that 
in this case the shot noise is not sufficiently removed. 
Actually, the shot-noise power spectrum amplitude exceeds 
the convergence power spectrum on scales below ~ 1 ar- 
cmin. In the following, therefore, we will only consider 
measured correlation functions on scales larger than 1 ar- 
cmin. The corresponding outputs will be labelled with 
small-scale smoothing below. 

Having produced surface density fields on all lens planes, 
1024 2 rays are traced backwards from the observer's point us- 
ing the multiple-lens plane algorithm (e.g. Schneider, Ehlers & 
Falco 1992). The initial ray directions are set on 1024 2 grids 
with a grid size of 0.25 arcmin, thus the total area covered by 
rays is 4.27 2 square degrees. We produced 36 realizations of 
the underlying density field by randomly shifting the simula- 
tion boxes in the direction perpendicular to the line-of-sight 
using the periodic boundary conditions of the A-body boxes. 
Note that the lens planes coming from the same box are shifted 
in the same way in order to maintain the clustering of matter in 
the box. 

We point out that second and higher-order statistics of 
point-source magnifications are generally ill-defined in pres- 
ence of caustic curves because the differential magnification 
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probability distribution asymptotically decreases as /i~ 3 for 
large /i (see Fig. ||). This is a generic feature of magnifica- 
tion near caustics and is thus independent of the lens model. 
Strong lensing effects on point sources near caustic curves give 
rise to rare, but arbitrarily high magnification values in the 
simulations, and therefore the variance of the measured statis- 
tics of fi cannot be defined. However, the smoothing proce- 
dure introduced above allows this problem to be removed be- 
cause it smoothes out high density regions in the dark mat- 
ter distribution and thus the fractional area of high magnifica- 
tion decreases. In reality, infinite magnifications do not occur, 
for two reasons. First, each astrophysical source is extended 
and its magnification (given the surface brightness-weighted 
point-source magnification across its solid angle) remains fi- 
nite. Second, even point sources would be magnified by a fi- 
nite value since for them, the geometrical-optics approximation 
fails near critical curves and a wave-optics description leads to 
a finite magnification (Schneider et al. 1992, Chap. 7). 




1 10 



Figure 2. Probability distribution of the magnification for 
our small- and large-scale smoothing simulations, assuming 
sources at redshift unity. The power law tail behaviour (oc fi~ 3 ) 
found in the small-scale smoothing indicates the existence of 
caustics, while for large-scale smoothing, no prominent tail is 
shown which suggests caustics do not play a noticeable role. 



3.2. Filtering 

The computation of correlation functions from numerical sim- 
ulations is mainly affected by two effects; on large scales 
by the finite box size of the dark matter simulation, and on 
small scales by the grid size used for computing the surface 
density field from the particle distribution. These boundaries 
set the limits for the validity of correlation functions mea- 
sured in numerical simulations. In other words, this means 



that measuring a correlation function on a given scale is rel- 
evant only if this scale falls within the range of scales de- 
fined by the simulation. As shown in the previous section, our 
method for computing the cross-correlation between n and k 2 
consists of first computing a three-point correlation function 
(k((P)k{4> + 6i)k(4> + 62)), and then identifying two of its 
three points. In such a case, one of the correlation lengths of 
the triple correlator becomes zero, thus necessarily smaller than 
the smallest relevant scales in any simulation. This prevents us 
from using any numerical simulation for directly comparing the 
results. 




Redshift 

Figure 3. Smoothing angle of the simulation as a function of 
redshift for the two ray-tracing schemes. In order to show the 
relevant quantities leading to the effective smoothing angle, we 
overplot the weighting function W(z) = D+(z) p^(z) (see 
Eqs. g] and 0]). 



In order to avoid this problem, and for comparing our an- 
alytical with numerical results, we will introduce an effec- 
tive smoothing into the theoretical calculations, such that each 
value of k at a given position 9 is evaluated by averaging the 
K-values in a disk of radius 6>g centred on 9. Indeed, the limit 
imposed by the grid size of the simulation gives rise to an un- 
avoidable smoothing-like effect which cancels all information 
coming from scales smaller than a corresponding smoothing 
scale 9s- For this purpose, we introduce a smoothed three-point 
correlator, 

z S mooth(02 - 0, 9 3 - 9A = ( k (9 1 )k{9 2 ) k {9 3 )) 9s (18) 

= J d9[j d9' 2 j do' 3 («(ei)«(ei)K(^)> 
x - OJWM - e 2 )w e M - 03) , 

where the function We s (9') is a normalised top-hat window of 
radius 9s- Introducing this smoothing scheme into the expres- 
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sion for (k(4>)k 2 (4> + 9)) yields 



3.3. Comparing (kk) and (kk 2 



Zsmooth(PSj 



dw- 



pI(w) f d 2 si f d 2 s 2 



■fft(w)J (2n) 2 J (2n) 2 
x l(s 1 9s)l{s 2 e s )l{\s 1 +s 2 \9 s ) 
si s 2 -«i - s 2 



x B s 



/fcM'/fcM' f k {w) 



where I(x) 



_ 2 Jife) 



Similarly, introducing the smoothing 



scheme into the two-point correlation function gives 

W smO oth(|02 - = ( K (^l) K ( 2))s s 



dw 



P 2 n(w) f d 2 s 
f 2 K (w) J (2n) 2 

|2 „is0 



P 



x |I(s6»s)|^e ls " . (20) 
The effective smoothing scale depends on two parameters: 

- the evolution of the apparent grid size of the simulation as 
a function of redshift, and 

- the radial selection function of the dark-matter field whose 
correlation function has to be measured. 

These quantities are plotted in Fig. ||. In order to use a 
unique smoothing scale valid on the final convergence map, 
we define the effective angular smoothing scale by 



dzw(z)e eiid (z) 



(21) 



where W(z) is the relevant normalised selection function along 
the line-of-sight. Measuring w s means probing the power spec- 
trum along the line-of-sight, weighted by p^iz). Therefore, we 
will use W{z) = D 2 + (z)p 2 K (z), where D + (z) is the growth 
factor. In a similar way, we will use W(z) — Di_(z)p\(z) for 
measuring z S mooth- The numerical values of the corresponding 
effective angles are presented in Table [[} 



small-scale smoothing 


large-scale smoothing 


^smooth #s = 0.40 


0s = 0.80 


Zsmooth 0s = 0.39 


9s = 0.78 



Table 1. Effective smoothing angles in arc minutes for w s and 
z s computed from Eq. [2l] as a function of simulation resolution. 



With the help of the filtering schemes introduced in the previ- 
ous section, we can now compare our theoretical predictions 
with correlation functions measured from the numerical simu- 
lations. We first compare the amplitude and angular variation of 
the two first terms of the Taylor expansion of the magnification 
separately. In the next section, we will then compare their sum 
to the total magnification fully computed from the simulation. 




8 (arcmin) 

Figure 4. Comparison between theoretical predictions and 
measurements from numerical simulations assuming sources 
at redshift unity. The upper and lower curves show (kk) and 
(kk 2 ), respectively. The points are measurements from the 
large-scale smoothing simulations, with the error bars show- 
ing the variance among 36 different realisations. The dotted 
lines show the analytical computations taking into account the 
smoothing scale of the simulation. The solid lines additionally 
include a cut in the power spectrum for cancelling the wave- 
lengths not covered by the simulation. The dashed line presents 
the same statistics without any smoothing. Obviously, the 
smoothing effects are crucial for the (kk 2 ) cross-correlation. 



The second important difference between analytical calcu- 
lations and measurements in numerical simulations is the finite 
box size effect. Indeed, the analytical correlation functions pre- 
sented above were computed taking into account all modes in 
the power spectrum. However, the finite size of the box used in 
the simulation introduces an artificial cutoff in the power spec- 
trum since wavelengths larger than the box size are not sampled 
by the simulation. This effect can also be taken into account in 
the analytical calculations by simply cancelling all the power 
on wavelengths with wave number k < fc mm - The boxes we 
use have a comoving size of 480 h Mpc which corresponds 
to fcmin = 0.013 /iMpc" 1 . 



In Fig. we overplot analytical and numerical results. The 
upper curve shows the autocorrelation function of k as a func- 
tion of angular scale. We plot in circles the average measure- 
ment from 36 realisations of the simulation, and the corre- 
sponding 1-er error bars to show the accuracy of the numerical 
results as a function of angular scale. The solid line shows the 
analytical prediction, including effective smoothing and an ar- 
tificial cut of the power at scales below fc m i n - The agreement is 
good on all scales. For comparison, the dotted line shows the 
result if we do not impose the large-wavelength cut, and the 
dashed line is the result if no cut and no smoothing are applied. 
In both cases, the deviations from the fully filtered calculation 
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Figure 5. The upper panel shows a comparison between the magnification autocorrelation measured from the simulation com- 
pared to the analytical estimation for a source redshift of unity. The circles show averaged measurements from 36 realisations of 
the simulation and the corresponding 1-a error bars. The dashed line shows the analytical estimation using the approximation 
\i fa 1 + 2k. The solid line shows the improvement given by the second-order term of the Taylor expansion of fi. The lower panel 
shows the relative contribution R^, both measured from the simulation and estimated using expansion terms of /i up to second 
order. In the right panel, each angular point shows three different measurements, taking into account the complete magnification 
maps (solid-line error bars), or maps where pixels were masked where the magnification value exceeds 8 or 4 (dashed-line error 
bars, from left to right). 



remain small since we are probing angular scales within the 
range allowed by the simulation. 



The lower curves in Fig. [4| show a quantity proportional to 
the second-order correction of the Taylor expansion, namely 
the correlation function (k k 2 ). In the same way as before, the 
circles show average measurements from 36 realisations, and 
the error bars denote the corresponding l-rr deviation. The pre- 
diction including smoothing and small-wavelength cut (solid 
line) shows a relatively good agreement given the expected 
accuracy of the bispectrum fitting formula, which is approxi- 
mately 15% (Scoccimarro & Couchman 2000). This time, in- 
cluding smoothing changes the amplitude dramatically, and 
this effect affects all scales (see the dashed line). As discussed 
before, this is expected since we are measuring a three-point 
correlator on triangles which have one side length smaller than 
the angular grid size of the simulation. Finally, as shown by 
the difference between the dotted and solid lines, cancelling the 
power on scales where k < k m i n again improves the agreement 
on large scales. 



The agreement between our analytical and numerical com- 
putations of (kk) and (k k 2 ) demonstrates the validity of the 
formalism introduced in Sect. Has well as the choice of the 
effective smoothing scale (Eq. [El]]) for describing the second- 
order term in the Taylor expansion of the magnification. 



3.4. Deviations from simulated ^-statistics 

We now want to investigate how well the second-order expan- 
sion describes the full magnification expression (§) which can 
be computed using maps of k, 7 and uo (a net rotation term 
which arises from lens-lens coupling and the lensing deflection 
of the light ray path; see Van Waerbeke et al. 2001b) obtained 
from the simulations (see Hamana et al. 2000 for more detail). 

Before doing so, we recall that the amplitude of the magni- 
fication autocorrelation measured from the simulation depends 
on the smoothing scale, as seen in Sect. since fi is nonlinear 
in the density field. Therefore, all the following comparisons 
are valid for a given effective smoothing length only. 

We further emphasise that two problems will complicate 
this comparison. First, our analytical treatment is valid in the 
weak-lensing regime only, i.e. as long as convergence and shear 
are small compared to unity, k <C 1, I7I <C 1. While most 
light rays traced through the numerical simulations are indeed 
weakly lensed, a non-negligible fraction of them will experi- 
ence magnifications well above two, say. Such events are re- 
stricted to small areas with high overdensities and thus af- 
fect the magnification statistics only at small angular scales. 
Second, a separate problem sets in if and where caustics are 
formed. The magnification of light rays going through caustics 
is infinite, and the magnification probability distribution near 
caustics drops like /i~ 3 for [i — > 00. As noted above, second- or 
higher-order statistics of /1 then become meaningless because 
they diverge. 
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Departures of the numerical from the analytical results will 
thus have two distinct reasons, viz. the occurrence of non-weak 
magnifications which causes the analytical to underestimate the 
numerical results on small angular scales; and the formation of 
caustics, which causes second-order magnification statistics to 
break down entirely. Both effects will be demonstrated below. 
They can be controlled or suppressed in numerical simulations 
by smoothing, which makes lensing weaker, or by masking 
highly magnified light rays or regions containing caustics. 

In Fig. ||, we plot with circles the autocorrelation func- 
tion (fJ,(4>) n{4> + ff)) measured from the large- and small- 
scale smoothing simulations in the left and right panels, respec- 
tively. The presence of caustics is more pronounced in the case 
of small-scale smoothing than in the large-scale smoothing 
simulations. The dotted line shows the theoretical prediction 
given by the first-order term of the Taylor expansion, namely 
4 (k(<P)k(<P + 9)). This yields a low estimate of the correlation, 
with a discrepancy of order 10% on large scales, and more than 
20% below a few arc minutes. 

As expected from the preceding discussion, this level of 
discrepancy also depends on the effective smoothing scale and 
can increase if simulations with a smaller grid size are used. 
Estimating the contribution of the two lowest-order terms of 
/i" -1 , we computed in Sect. 2.3 a lower bound to this discrep- 
ancy for a real case without smoothing, and found it to reach 
a level of 25% at large scales, and above 30% below a few 
arc minutes. The smoothed results taking the additional con- 
tribution of the second-order term into account are plotted as 
solid lines, and give a much better agreement, as expected. To 
quantify this in more detail, the lower panels of the figure show 
several contributions compared to the first-order term, i.e. to 
4(kk). 

- The symbols show the additional amplitude of the magni- 
fication statistics measured from the simulation, compared 
to the first-order term also obtained from the simulation, 



in 

= 4( K (0M0 + 0)) r 



- 1 



(22) 



The error bars indicate the 1-cr deviation across 36 realisa- 
tions. 

The solid line shows the contribution of the second-order 
relative to the first-order term computed from the analytical 
expression including the effective smoothing, 



,„„i, _ 2a(K(cj))K 2 ((l) + 0)) es 



R 



with a = 2. 



(k(0)k(0 + 0)>0 s 



(23) 



In each case, we use the appropriate reference for 
(n(<p)n{<p + 9)), i.e. the numerical measurement in the first 
and the analytical estimation in the second case. Indeed, the 
measurement of (k k) from the simulation agrees with the ana- 
lytical estimation within some uncertainty, which is due to nu- 
merical effects like the finite number of dark-matter boxes used 
for simulating the light cone. It introduces a bias into our com- 
parisons which is impossible to separate from the real offset 
due to all higher-order terms of the Taylor expansion that were 



not taken into account. The two contributions plotted in Fig. || 
are thus of different nature, but are suitable for a relative com- 
parison. 

As the lower panel of the large-scale smoothing simula- 
tion shows, the simple 4 (k k) estimate of the magnification 
misses 20% of the real amplitude near one arc minute. This 
discrepancy almost vanishes after adding the contribution of 
the second-order term, which gives at all scales a final agree- 
ment on the per cent level: the additional amplitude reaches 
19% at the smallest scales of the figure, compared to a value 
of 20% given by the simulation, and agrees within better than 
one per cent on larger scales. Therefore, taking into account the 
2 a (kk 2 ) correction allows the accuracy to be increased by a 
factor of ~ 20 compared to the approximation 4(kk), in the 
case of our large-scale smoothing simulation. On the largest 
scales, between 6 and 30 arc minutes, the agreement even im- 
proves. Above these scales, the numerical results do not allow 
any relevant comparison because the number of available in- 
dependent samplings corresponding to a given separation de- 
creases. On scales below a few arc minutes, the offset between 
the measured points and the analytical estimate gives the ampli- 
tude of all higher-order terms neglected in the Taylor expansion 
of the magnification. As we can see, their contribution is on the 
one per cent level for the large-scale smoothing simulation. 

The curves shown in the right panel demonstrate how the 
use of a smaller smoothing scale increases the discrepancy be- 
tween the analytical and the numerical results. The fraction of 
non- weakly magnified light rays increases, and caustics appear 
which give rise to a power-law tail in the magnification proba- 
bility distribution. We investigate the impact of the rare highly 
magnified light rays by masking pixels where the simulated 
magnification exceeds 4 or 8, and show that caustics have no 
noticeable effect on the amplitude of the magnification autocor- 
relation function determined from these simulated data. Note, 
however, that the impact of the caustics depends on the source 
redshift. The higher the redshift, the more caustics appear, and 
the larger is their impact on the correlation amplitude. 

Imposing lower masking thresholds removes a significant 
fraction of the area covered by the simulation, changing the 
spatial magnification pattern and thus the magnification auto- 
correlation function. The corresponding measurements are rep- 
resented by the dashed error bars in the lower right panel of 
Fig. |[ We note that the error bars of i?™ m computed with 
the small-scale smoothing simulation become larger at small 
scales compared to the lower left panel. This reflects the fact 
that second-order magnification statistics are ill-defined once 
caustics appear. In the next section, we will investigate similar 
smoothing effects on cross-correlations between magnification 
and dark matter fluctuations. These quantities are not affected 
by problems of poor definition when the smoothing scale be- 
comes small, and therefore do not show larger error bars at 
small scales when the smoothing scale decreases. 

These comparisons show that the approximation /j, w 1+2k 
misses a non-negligible part of the total amplitude of weak- 
lensing magnification statistics. The formalism introduced in 
Sect. U allows second-order corrections to be described with or 
without smoothing of the density field. This provides a better 
description of the correlation functions, but still gives a lower 
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amplitude than the simulation results. As we noticed, the ana- 
lytic computation based on the Taylor expansion is sufficiently 
accurate only in the weak lensing regime. In reality, however, 
the strong lensing, which can not be taken into account in the 
analytic formalism, has a significant impact on the magnifica- 
tion correlation especially at small scales as shown in the small- 
scale smoothing simulation. Therefore, one should carefully 
take the strong lensing effect into consideration when one inter- 
prets the magnification related correlation functions. However, 
we will see in the next section that counts-in-cells estimators 
are less affected by the strong lensing than correlation func- 
tions and thus enable better comparisons of observations with 
results from simulations. 

4. Applications to quasar-galaxy correlations 

As a direct application of the formalism introduced previ- 
ously, we now investigate the effects of second-order terms 
on a well-known magnification-induced correlation, namely 
the quasar-galaxy cross-correlation (the results can also be ap- 
plied to galaxy-galaxy correlations induced by magnification; 
Moessner & Jain 1998). In order to estimate cosmological pa- 
rameters from this kind of correlations, we then suggest the use 
of a more suitable estimator using counts-in-cells rather than 
two-point correlation functions. It has the advantage of making 
the observational results more easily reconciled with the ones 
from numerical simulations. 

4. 1. Formalism and correcting terms 

The magnification bias of large-scale structures, combined 
with galaxy biasing, leads to a cross-correlation of dis- 
tant quasars with foreground galaxies. The existence of this 
cross-correlation has firmly been established (e.g. Bemtez & 
Martinez-Gonzalez 1995; Williams & Irwin 1998; Norman & 
Impey 1999; Norman & Williams 2000; Bemtez et al. 2001; 
Norman & Impey 2001). Menard & Bartelmann (2002) showed 
that the Sloan Digital Sky Survey (York et al. 2000) will allow 
this correlation function to be measured with a high accuracy. 
Its amplitude and angular shape contain information on cosmo- 
logical parameters and the galaxy bias factor. Thus, it is impor- 
tant to accurately describe these magnification-related statistics 
in order to avoid a biased estimation of cosmological parame- 
ters as well as the amplitude of the galaxy bias. 

As shown in Bartelmann (1995), the lensing-induced cross- 
correlation function between quasars and galaxies can be writ- 
ten as 

w QG (6) = (d Q soW>)<5 ga iW> + 0)) 

= (^ Q - 1 (0),5 gal (0 + 0)} . (24) 

Using the above formalism, we can expand the effective mag- 
nification fluctuation 5fi a ~ 1 up to second order and find the 
correcting term: 

w QG (0) = 2 (a -I) [( K 6 sal ) + a ( K 2 S gal )] . (25) 

The second term is proportional to a (contrary to the factor 2a 
in Eq. |Q|), since there is only one contribution of the magni- 
fication. Therefore, the expected effects will be roughly half 



of those on the autocorrelation of the effective magnification 
seen in the previous section. Assuming a linear bias b between 
galaxies and dark matter, the cross-correlation between (5 ga i and 
K 2 can be written as 



<<5 gal (0) K 2 (0 + 0)> 
pl(w)p s (w) 



dw 



f 4 K H 



x B s 



•*>'2 



b(s DM (<P) K 2 {cf> + e)) 

d 2 Sl f d 2 s 2 

(27T) 2 

-Si - S 2 



(27T) 2 



(26) 



Jk(w) fk(w) f k (w) 

where ps{w) is the normalised distance distribution of the 
galaxies. For this example, we will use 



p s (z) dz 



(3z 2 



; iW) 



exp 



dz 



(27) 



with (3 = 1.5 and zq = 0.3. 

The results are shown in Fig. || As we can see, previous 
estimates using the approximation fi « 2 k missed approxi- 
mately 15% of the amplitude on small scales for quasars at 
redshift unity. Using quasars at redshift 2, these effects reach 
up to 25%. These offsets, which are only lower limits, would 
lead to biased estimates of SIq or b, for example. 
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Figure 6. The upper panel shows the amplitude of the nor- 
malised quasar-galaxy correlation WQG/2(a — 1) as a function 
of angular separation. We show the first two terms of the Taylor 
expansion of this correlation, namely {8(4>)k(4> + 6)) (dashed 
line) and 2(5(4>)k 2 (4> + 6)) (dotted line), using a source red- 
shift of unity. The sum of these two terms is shown as the solid 
line. The lower panel details the relative contribution Rg^ and 
of the second-order term for different source redshifts, namely 
z = 1, 1.5, 2 and 3 from bottom to top. 



As for the magnification autocorrelation, we can compare 
our theoretical estimates against numerical estimations. We can 
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first introduce a coefficient i?^ describing the accuracy of our 
second-order correction: 



R Sfl (0) 



{5{4>)^{4> + e)) 



(28) 



(5(4>)k(<I> + e)) ■ 

We plot the results in Fig. [7| Note that contrary to the magnifi- 
cation autocorrelation, this quantity does not suffer from poor 
definition, even without smoothing. The difference can be seen 
by the same size of the error bars between the two simulation 
results at small scales, whereas they were larger in the case of 
(fifj,) for the small-scale smoothing simulation (Fig. ||). The re- 
sults for Rs/j, are very similar those obtained for R^: for the 
large-scale smoothing ray-tracing we find very good agreement 
which reaches the one percent level on small scales. However, 
when the smoothing length decreases, we see from the small- 
scale smoothing outputs that we are missing a part of the total 
amplitude on small scales, which shows that higher-order terms 
play a non negligible role on those scales. 



20% 




05% - 



00% 
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Figure 7. Comparison of Rs^ given by the theoretical calcula- 
tion and measured from the numerical simulation. The results 
are shown for the large- and small-scale smoothing simulations 
from bottom to top. They show that a second-order descrip- 
tion of the (fj,S) cross-correlation gives good results when the 
smoothing is large, but misses some part of the total amplitude 
in the case of our small-scale smoothing ray-tracing. 



4.2. Using count-in-cells estimators 

For precisely estimating cosmological parameters as well as the 
amplitude of the galaxy bias, it is necessary to employ theoret- 
ical magnification statistics that closely describe the observ- 
ables. However, we have seen in Sect. |5] that analytical esti- 
mates as well as numerical simulations have intrinsic limita- 



tions and prevent us from accurately describing usual n-point 
correlation functions related to magnification statistics. 

Besides, it is possible to focus on another estimator closely 
related to correlation functions, namely a count-in-cells esti- 
mator, which naturally smoothes effects originating from the 
density field and can thus more easily be reconciled with nu- 
merical simulations. So far, quasar-galaxy or galaxy-galaxy 
correlations have been quantified measuring the excess of 
background-foreground pairs at a given angular separation. 
Instead, we can correlate the amplitude of the background and 
foreground fluctuations, both measured inside a given aperture. 
We will therefore introduce a count-in-cells estimator, 



Wqg(6) = (<5QSoO)<SgalO))fl 

( ( 5 M a - 1 (0) ( 5 gal (0)) 6 



(29) 



where the subscript 8 indicates averaging of <5qso(0) an d 
S ga i((f>) inside a cell of radius 8. In practice, this estimator is 
intended to be applied to galaxy-galaxy rather than to quasar- 
galaxy correlations, since the average angular separation be- 
tween bright distant quasars is of order one degree for cur- 
rent surveys, thus averaging the source counts inside cells with 
radii of several arc minutes will not be relevant. Using galaxies 
as background sources, this limitation occurs at much smaller 
scales. 

Using a first-order Taylor expansion for the magnification, 
the new estimator wqg(6) can be written 



2 (a- 1)6 



where I(x) 



dw 



p K (w)ps(w) f sds 
2tt 



X P 



/fc(w)' 



\l(sd)f 



(30) 



_ r, JlQ) 



This expression differs from the 2- 



point correlation function (|9J) by its Fourier-space filtering of 
the power spectrum. The additional smoothing wipes out the 
power on scales smaller than the physical scale corresponding 
to the angular smoothing scale 8. For any observational result 
to be compared to a numerical simulation, 8 and the smoothing 
scale used in the simulation will have to be carefully adapted 
to each other and to the redshift distribution of the foreground 
galaxy distribution. 

Note that gravitational lensing by the foreground galax- 
ies themselves is entirely irrelevant here. The angular scale on 
which galaxies act as efficient lenses is on the order of one 
arc second and below, much smaller than the angular scales 
we are concerned with. Moreover, the probability for a quasar 
to be strongly lensed by a galaxy is well below one per cent. 
Bartelmann & Schneider (1991) demonstrated this point ex- 
plicitly by including galaxies into their numerical simulations 
and showing they had no noticeable effect. 

5. Conclusion 

As surveys mapping the large-scale structure of the Universe 
become wider and deeper, measuring cosmological parameters 
as well as the galaxy bias with cosmic magnification will be- 
come increasingly efficient and reliable. Therefore, an accurate 
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theoretical quantification of magnification statistics becomes 
increasingly important. 

Previous estimates of cosmic magnification relied on the 
assumption that the magnification deviates sufficiently little 
from unity that it can be accurately approximated by its first- 
order Taylor expansion about unity, i.e. fj, m 1 + 2 k. In this pa- 
per, we have tested the validity of this assumption in the frame- 
work of magnification statistics, by investigating the second- 
order terms in the Taylor expansion of /i. We have shown that: 

- Second-order terms can be related to the cross-correlation 
between k and k 2 , 

- their importance increases as the number-count function of 
the background sources steepens, i.e. as a increases, 

- their amplitude is not negligible: for the magnification au- 
tocorrelation, their contribution is typically on the order of 
30%-50% at scales below one degree. Therefore, previous 
estimates of cosmic magnification were systematically low. 

For testing our theoretical calculations, we have compared 
our results to magnification statistics found in numerical sim- 
ulations by performing ray-tracing experiments in a very large 
iV-body simulation. We have first checked the validity of our 
formalism describing the correlation (kk 2 ), and demonstrated 
the importance of including an effective smoothing into the an- 
alytical calculations. Indeed, [i is nonlinear in the density field 
and the amplitude of magnification statistics measured from 
numerical simulations depends therefore on the available reso- 
lution. 

Using a simulation with an effective smoothing scale of 
0.8 arc minutes, we found that our second-order formalism is 
accurate to the percent level for describing magnification au- 
tocorrelations. Compared to previous estimates, this improves 
the accuracy by a factor of rj 20. For smaller effective smooth- 
ing scales, the contribution of third- and higher-order terms be- 
comes important on scales below a few arc minutes. 

Finally we have applied our formalism to observed correla- 
tions, like quasar-galaxy and galaxy-galaxy correlations due to 
lensing. We have shown that second-order corrections increase 
their amplitude by 15% to 25% on scales below one degree. 
These correlations are valuable tools to probe cosmological 
parameters as well as the galaxy bias. However, even includ- 
ing our correcting terms, analytical or numerical estimates of 
magnification statistics can only provide lower bounds to the 
real amplitude of the correlation functions in the weak-lensing 
regime. Thus, we propose using count-in-cells estimators rather 
than correlation functions since the intrinsic smoothing in de- 
termining counts-in-cells allows the observational results to be 
more directly related to those obtained in numerical simula- 
tions. 
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Appendix A: Bispectrum and Non-Linear 
evolution 

The bispectrum can be estimated using second-order perturba- 
tion theory. Indeed, an expansion of the density field to second 
nonlinear order as 



6(x) ra <5 (1) (x) + 5 {2 \x) , 



(A.l) 



where 5^ is of order (S^) 2 and represents departures from 
Gaussian behaviour, yields the bispectrum 



(616263) 



(6^6^>6^) + (6\ 



(l)jr(lM 2 ) 



cyclic terms (231, 312) 



(A.2) 



The first term in Eq. (A.2) vanishes because the density fluc- 
tuation field is Gaussian to first order, hence the third moment 



of is zero. Thus, the leading term in Eq. (ft..2|) is of the or 



der of (^ 1) 4 1) 4 2) ) and 

can be quantified using second-order 

perturbation theory. 

The bispectrum Bg (ki, k 2 ,k 3 ) is defined only for closed 
triangles formed by the wave vectors ki, k2, k 3 . It can be ex- 
pressed as a function of the second-order kernel F(ki, k2) and 
the power spectrum 



S 5 (ki,k 2) k 3 



2F(k 1) k 2 )P(fe 1 )P(fc 2 ) 
2 F(k 2) k 3 ) P(k 2 )P(k 3 ) 
2F(k 1) k 3 )P(fc 1 )P(fc 3 ) 



(A.3) 



For describing the bispectrum on all scales, we use the fitting 
formula derived by Scoccimarro & Couchman (2000) for the 
non-linear evolution of the bispectrum in numerical simula- 
tions of CDM models, extending previous work for scale-free 
initial conditions. In that case, we have 



F(ki,k 2 ) = 



5 
7 

Iki 

2 k x k 2 



a(n, fci)a(n, fc 2 ) 
k 2 ( ki fc 2 



with 

a(n, k) = 

b(n, k) = 

c(n, k) = 



1 + Oz 



ki • k 2 

k\k 2 



fci 

fc 2 
2 



b(n, ki)b(n, fc 2 ) 



c(n, ki)c(n, fc 2 



(A.4) 



\z) [0.7Q 3 (n)] 1/2 (g/4) 



n+3.5 



1 + (g/4)«+ 3 - 5 
1 + 0.4 (n + 3) q n+3 

1 + g«+3.5 

1 + 4.5/ [1.5 + (n + 3) 4 1 (2q) n + 3 



1 + (2 9 )"+ 3 - 5 



(A.5) 



and q = k/kn^z), where 47rfc^ L PL(fcNL) = 1, and Ph(k) is 
the linear power spectrum at the desired redshift. The effective 
spectral index is taken from the linear power spectrum as well. 
The function Q 3 (n) is given by 



Qa(n) 



(4-2™) 



(1 + 2™+!)' 

For more detail, see Scoccimarro & Couchman (2000). 



(A.6) 
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